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Friedberg, Lee and Zhao proposed a method for effectively evalu- 
ating the eigenenergies and eigen wavefunctions of quantum sys- 
tems. In this work, we study several special cases to investigate 
applicability of the method. Concretely, we calculate the ground- 
state eigenenergy of the Hellmann potential as well as the Cornell 
potential, and also evaluate the energies of the systems where 
linear term is added to the Coulomb and harmonic oscillator po- 
tentials as a perturbation. The results obtained in this method 
have a surprising agreement with the traditional method or the 
numerical results. Since the results in this method have obvious 
analyticity compared to that in other methods, and because of the 
simplicity for calculations this method can be applied to solving 
the Schrodinger equation and provides us better understanding of 
the physical essence of the concerned systems. But meanwhile ap- 
plications of the FLZ method are restricted at present, especially 
for certain potential forms, but due to its obvious advantages, it 
should be further developed. 



I. Introduction 

Recently, Friedberg, Lee and Zhao (FLZ) have proposed a method Q which is very powerful 
and useful for solving the Schrodinger equation, especially for the potentials which have strong 
couplings. 

In their works [^, applicability and characteristics of the method were discussed, on both 
the physical and mathematical aspects. 
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As well known, only very few potential forms can lead to analytic solutions of the Schrodinger 
equation, whereas for the majority of potentials which are of phenomenological significance in 
real physical world, one can only expect numerical results. It reduces the ability for further 
analyzing the problem and getting insight into the physics essence. Therefore an almost-analytic 
solution would be very welcome for studying the physics behavior of the system. 

In quantum mechanics, the WKB approximation is an important method to study the 
potential-barrier tunneling and quantization conditions for bound states. In fact, the WKB 
approximation is an semi-classical approximation where one can write the wavefunction -0 as 
■ip{x) = exp{^a{x)) and expand a{x) according to orders of h. Generally, only the first two 
terms are retained It is confirmed to be useful for tunneling rate evaluation, but for the 
bound states, it can only offer a rough estimation. 

Comparing with the WKB approximation, FLZ proposed that the wavefunction can be 
written as 

<A(g) = e-^(«), (1) 

where g is a set of any coordinates. Obviously, the WKB wavefunction exp{^a{x)) is an oscil- 
lating form which corresponds to real particles, so that more applicable to the tunneling case, 
whereas, the FLZ wavefunction is a damping form and more suitable for dealing with bound 
states. 

In this work, we would like to investigate the applicability of the FLZ method in terms of a 
few potentials which have significance for practical physics problems. We first, in terms of the 
FLZ method, evaluate the eigenenergy of the ground state for the Cornell potential ^ and the 
Hellmann potential [^, and secondly we study applications of this method to the perturbation 
calculations. Concretely, we consider a linear term as a perturbation to the hamiltonians which 
contain either the Coulomb potential(the Cornell potential) or harmonic oscillator potential, 
because the first one has importance for spectra of heavy quarkonia and the second has a 
rigorous analytic solution. 

We also employ the Dalgarno-Lewis technique to obtain the corresponding perturbation 
results up to the third order (for the Cornell potential). Comparing the results with that 
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obtained in the FLZ method, one can find that they coincide with each other order by order. 
However, obviously the procedure of derivation in the FLZ method is much simpler. 

It is noticed that the FLZ method cannot be applied to the cases for linear potential aiq, or 
the potential (not for a perturbation). We will briefly discuss this issue in the last section 
and also concern other restrictions to applications of the method. 

This paper is organized as following, after this introduction, we briefly introduce the con- 
cerned aspects of the FLZ method and in Sec. Ill, we discuss the Hellmann potential case and the 
Cornell potential case, then in Sec. IV, we apply this method to the perturbation calculations. 
The last section would be devoted to the conclusion and discussion. 

II. Brief introduction to the FLZ method 

For a completeness and convenience of readers, we briefly introduce the Friedberg-Lee-Zhao 
method in this section. 

The Schrodinger equation can be written as 

Hcl){q) = E4>{q), (2) 

where the hamiltonian is 

H = ^V^ + V{q) and = E S 

i=l 

and here for simplicity we take the effective mass /i to be unity. Friedberg, Lee and Zhao suggest 
to write the wavefunction of the ground state in the form of eq.(|l]) and moreover, they set 

V{q)=gMQ), (3) 

where is a scale factor and corresponds to the strength of the effective coupling, for example in 
the QED Coulomb case, it is Ze^ and for the QCD Coulomb case, it is Oe// |||]- The advantage 
of pulling out the coupling parameter will be shown in our later examples. 

Choosing as an expansion parameter, one can expect from the expression of eq.(|l|) that 
the larger g^ is, the faster the expansions of s{q) and E converge. Thus one can have expansions 
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as 

E = g'Eo + g'-'Ei + g'-^'E2 + ..., (4) 
s{q) = g^'so + g'^-hi+g"'-^^S2 + ..., (5) 

where /, m, i.j are positive integers and to be determined according to a rule that under hmit 
of strong interaction, 

n + 2 ^ ' 

where n is the power of the leading term in the potential, for example, n = 2, I = k/2 for the 
harmonic oscillator and n = — 1, I = 2k for the Coulomb potential. To avoid difficulties in the 
WKB approximations such as the matching conditions at the turning points, it is required that 
£"0 and V do not appear in one equation and other equations are solvable. Thus one has 

I < k : 2m = k, 2m — j = I, i = j, (7) 

I > k : 2m = I, 2m-j = k, i = j, (8) 

where k is any positive integers which can make / to be an integer. 
For example, in the simplest case, l = m = i = j = l, one has 

{Vsof = 2v, 
Vso-Vsi = ^V^so-^o, 
Vso-Vs2 = ^{V^si-{Vsif)-Ei, 

(9) 



It is noted that in the first equation, Eq does not appear, and only v exists, whereas in the 
second, only £^0) but v is absent. So in this method, it is not needed to distinguish between 
E < V or E > V cases. 

From the first equation of @, one can obtain Vsq and then substitutes it into the second 
equation. Here it is supposed that the potential is a function of r. For the S-states, = 



^ + faf'thus 



d 

Sm\r—>0 



dr 



must be zero to insure (V'^Sm+i)|r^o finite, where m is the order index. With this condition, we 
obtain Eq and Vsq simultaneously from the second equation. Then we will continuously achieve 
Em-i and Vsm from the m-th equation. 

With the boundary condition (/)(0) = 1 or s(0) = 0, we also obtain the wavefunction up to a 
normalization factor. 

Later we will show that this method is very powerful for evaluating the eigenenergy of the 
ground state. 



III. Applications to the Hellmann potential and the Cornell potential 

(1) The Hellmann potential is a superposition of a Coulomb piece and a Yukawa piece 

V{r)=g\— + -e-^'), (10) 
r r 

which cannot be analytically solved and here we deliberately pull out a coupling constant . 
In physics, the Coulomb potential usually originates from a single massless photon or gluon 
exchange, but the Yukawa potential is due to a massive particle exchange and can be induced 
by a form factor at the effective vertices. Such a potentail may play an important role in 
phenomenology. Now let us evaluate the ground energy in terms of the FLZ method. We have 
the expansions for s and E similar to (^, and here / = 4, i = j = m = 2, then the set of 
equations is recast as 

[Vsof = -2Eo, 

Vso-Vsi = Ivho - - + -e-^'- - Eu 
2 J. J. 

Vso-Vs2 = ^(Vsi)2 - iv^si -^2, 

VS0-VS3 = -Vsi- Vs2 + ^V^S2-^3, 

(11) 



The general expression for the even order of n > 1 is 

n-l ^ ^ 

Vso • Vs2n = - X! ^^rn ' V S2n-m " -^[^ Snf + -V^S2n-l " -E'2n, (12) 



2' 

m=l 



and for odd order, it is 



1 

Vso • Vs2n+1 = - X! ■ '^S2n+l-m + ^^^"^8211 " -E'2n+1- (13) 

m=l 



Thus we obtain 



So = \/-2Eqj 



where we only keep the the positive root. The second equation is written as 

dr r r r r 

It is required that the left side of the equation and Ei must be finite at r = 0, thus a natural 
condition 

^/-2Eq -A + Be-^'' ^0, as r ^ 0, 
is enforced. One obtains immediately 

EQ = ^{A-Bf and SQ = \A-B\r. (15) 

As a matter of fact, the key point of the FLZ method is to obtain the first-order eigenenergy by 
removing the singularity in the equation and it is similar to the secular equation for obtaining 
eigenenergies in quantum mechanics. 

Because all V'^si (i=l,2...) are finite, so if we expand 

(j) {€) (i) 

(X j2 ^ n-Ul ^ 



n+l I , "-1 I {i) I (i) I (i) 2 



all and ai must vanish, thus 

— > 0, as r — > 0, (m > 1). 

dr 

But So is an exception which is linearly proportional to r, indeed, V'^sq is singular at r — > 0. It 
is exactly the condition to determine the value of so as done in eq. (|l^) . 

As we go on with the same strategy we obtain other E[s and one can find that all superficial 
singularities in the expressions are cancelled out automatically. We write the eigen-energy of 
the ground state as 

^ = - ^) V - BCg^ + + 0(-L), (n > 1). (16) 



Obviously as the coupling is sufficiently large, the higher orders can be safely neglected. 
Moreover, one can notice that as A = B \i blows up and we will discuss this issue later. 

To make sense, we take a set of the parameters to testify the results. We can numerically 
solve the Schodinger equation and obtain the numbers. Of course, we suppose that the numerical 

results arc precise. A comparison of the numerical results with that with the FLZ method is 
shown in Table 1. 



g 


Numerical result 


FLZ 


1 


-1.0848 


-0.75 


2 


-11.3687 


-11.25 


3 


-48.805 


-48.750 


5 


-336.757 


-336.750 



Table 1. The comparison of the numerical results with that in terms of the FLZ 
method where only the first three terms in the expansion are taken. The parameters 
are set as A = 2, B = C = 1. 

It is noted that as 5 = 3 ~ 5, the sum of the first three terms from the result in the FLZ 
method are perfectly consistent with the numerical results. 
(2) The Cornell potential. 
The hamiltonian is 

H = - M + Kr, (17) 

for using the FLZ method let us transform it into another convenient form 

H = ^V^ + g\— + n'r), (18) 

where g"^ is pulled out. The set of equations is similar to that for the Hellmann potentail as 

^ — ^e~^^ is replaced by ^ -|- nr. 

{Vs^f = -2Eo, 

1 9 1 , 
Vso-Vsi = -V^so-- + K'r-Ei, 

2 r 

Vso-Vs2 = ^{Vsif -^Vhi-E2, 

1 9 

(19) 
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The following procedures are exactly the same as in last subsection. We obtain 




(20) 



One can notice that in the hamiltonian ([l^), we set the mass /i = 1 for simplicity. When we put 



back the mass, the Schrodinger equation is 



[— V2-^ + Kr]|V'>= E\^> 



(21) 



the equation can be re-written as 



2 r 



(22) 



Therefore, the expansion parameter is replaced by ng^. Even though in the practical case, 
the effective coupling Oe// is about 0.4~0.5 being small, for heavy quakonia such as T, the 
reduced mass fi ~ 2.5 GeV, and fxg^ is not too small, the method can give result which is close 
to the numerical result. Our numerical calculation indeed shows that for the bb quakonium T, 
the FLZ method works very well, but for the cc system J ftp where /i ~ 0.75 GeV, the method 
fails to give a reasonable solution. 

Now let us turn to the perturbation cases. 

IV. Perturbation in the FLZ method 

(1) First we discuss a very simple case where a linear term is attached to the harmonic 
oscillator perturbation, the hamiltonian is 



for simplicity, here we still take the mass /i = 1. The accurate ground state energy is E = 
hd ~ Thus we would use this result to testify the applications of the FLZ method for 

dealing with a perturbation and in this situation a is the expansion parameter. According to 




(23) 



the general rule and eq.(P), we have 
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(P' dso dsi 

-r~2^o - 2-^-3— = 2Eo, 
ax^ dx dx 

d^si Jsi2 ^dsods2 

ax^ dx dx dx 

d^S2 _ ^dso_ds3_ _ ^ds^ds2 _ 

dx"^ dx dx dx dx 

(24) 

From the first and the second equations, 

^ ds\ 
dx 

setting X = 0, we obtain Eq = ^ and ^ = 0. As we go on calculating, we can achieve 

-1 9 

El = 0, E2 = 0, Ez = — a^ and Ei = E^ = ... = 0. 
So = -x^ , Si = 0, S2 = ax, ... 



So finally, 



^1 la2 

E = -g ^. 

2^ 252 

This energy is precisely the same as the exact solution. 



(2) The Coulomb potential with a linear term as the perturbation. 

In this subsection, we treat the linear term Kr as a perturbation to the Coulomb piece. The 
hamiltonian is 

-1 

H = — - — + Kr. (25) 
2 r 

The difference between eq.(p5|) with eq.(^) is the position of the coupling g^. Thus we have a 
set of equations as 

(Vsof = -2Eo, 

Vso-Vsi = lv^so---Ei, 
2 r 
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Vso-Vs2 = -^{Vsif-^Vhi + Kr-E2, 
Vso-Vs3 = -Vsi-Vs2 + ^Vh2- Es, 



(26) 



again, the difference of the set from the set (p!9|) is the position of nr which appears in the third 
equation, but in the set ([l9|), it exists in the second equation. Solving the equations according 
the normal procedure, we have 

1 3 —3 2 27 

-^0 El = E2 = 0, £"3 = -K, E4 = E^ = 0, Eq = —K , Ey = Es = 0, Eg = —k , ... 

si =r, si = 0, S2 = \Kr^, S3 = 0, 54 = 7r4, S5 = TT^r^ , ■■■ 
Thus the energy of the ground state is 

E = ^9' + l^g~' - l^-^g-' + f ^V^^ + 0(-V^°) + .... (27) 

If we set K = k! , we recover the expression ([20|). The first term ( ^7] ) is the energy caused by 
the Coulomb potential alone. It is amazing to notice that in this subsection, we treat kt as a 
perturbation to the Coulomb potential, while in section 111.(2), we deal with the Coulomb and 
linear pieces as a whole potential, then the two results are the same to the order 0(k^(7~^^). 

For a better comparison, let us calculate the contribution from the linear term order by order 
in the traditional perturbation theory. The first order perturbation result can be easily obtained 
in the traditional method as 



El =< i?i,o|Kr|i?i,o >= 77^, (28) 



3k 
2? 

where Ri^q = 2g^e~^^^. 

For evaluating the second order and even higher orders of corrections to the eigen-energy, 
the calculation procedure in the method given by most textbooks of quantum mechanics is too 
complicated to carry out practical calculations. It fails because in the traditional method, all 
energy states including the continuous and discrete spectra must be summed over. Dalgarno 
and Lewis introduced a technique |0, which allows us to calculate the perturbation contribu- 
tions without carrying out the tedious summation and obtain results at any given order. The 
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technique, in fact, is based on the hypervirial theorem |8| and ah details about the technique 
can be found in ref.|^. 

By this technique, we can immediately obtain the second and the third order corrections to 
the energies of the ground state of the Coulomb potential as 

_ -3a^2 27 

and this is exactly that we derived in the FLZ method. 
V. Conclusion and discussion 

Prom Sees. Ill and IV. where we discussed solutions for the Hellmann potential and the 
Cornell potential which has importance to the heavy quarkonium, and the perturbative cases in 
the FLZ method. First we find that for strong coupling g at the Hellmann potential the results 
obtained in terms of the FLZ method are very close to the numerical results, in fact the larger 
g is, the closer to the real solution the FLZ result is. We can conclude that the FLZ method 
applies perfectly for strong interactions. Moreover, this method allows us to have analytic forms 
for both eigenenergies and wavefunctions, so it has obvious advantages for physical analysis. 

Then we use the FLZ method to calculate the energy corrections due to a linear perturbation 
to the harmonic oscillator and find that the result is exactly the same as the precise solution. 

We employ the FLZ method to calculate the energy for the Cornell potential and pretend 
that the coupling is strong. Then alternatively, we treat the linear piece as a perturbation 
to the Coulomb piece. We calculate the correction to very high orders in the FLZ method 
and then we employ the Dalgarno-Lewis technique to repeat the calculation, amazingly, we 
find that all the three methods give the same results (precisely), even though they start from 
completely different points. This amazing consistence confirms the applicability of the FLZ 
method. Moreover, the calculations with the FLZ method are much simpler than with others 
and one can avoid the tedious and complicated integration and summation. The solution of 
both energy and wavefunction have almost-analyticity. In fact, from our pedagogical examples, 
one can believe that the FLZ method can apply to much more complicated perturbation terms, 
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such as the L-S couphng, tensor potentials and the relativistic corrections etc. 

On other side, at present, the FLZ method can only effectively apply to deal with the ground 
state. FLZ showed Q that for a special case the method can be applied to evaluate the eigen- 
energy of the first excited state, but not for general cases yet. How to generalize this method to 
calculate excited states in general is an interesting task and should be further investigated. 

The second serious drawback of the method is that the method cannot apply to the potential 
where only positive power terms exist and the leading order is the linear one. The reason can 
be understood as the following. If only the linear term exists, for example, all the equations in 
the (^) become trivial, so that one cannot determine Ei from them and the method fails. There 
could be some way to remedy this problem, and we will present some possible ways to solve in 
our later work. 

It is also shown in our expression for the Hellmann potential, the resultant energy is related to 
A — B, which sometimes appears at denominator. As A ^ B, the corresponding term diverges. 
Turning to the original expression of the Hellmann potential, as A = B, it just is the case of a 
linear term as the leading order. Because there exists an exact solution for the S-states in the 
case of the linear potential Kr and it is the well-known Airy function, therefore the superficial 
divergence is associated to the problem in the FLZ method. We will further investigate this 
issue in our future work. 

Moreover, it requires a strong coupling if one wants to have a more accurate result for a 
certain potential. For example, in Sec. HI (2), for the Cornell potential, we pretend or assume 
that the coupling is strong, i.e. g is large, however, in reality, the coupling in the Cornell 
potential is Oe// = ^f^ ~ 0.3 to 0.5 and not large at all. But as shown above, the real expansion 
parameter is fig^ instead of g"^, so as long as /i is large enough, fj,g^ is also large and the FLZ 
method applies. This is the case for the bb quarkonium T, but does not apply to the case of J/ip. 
Thus we, in general, so far cannot expect to simply calculate the charmonium ground states in 
this method yet. 

On other side, the method has obvious advantages as discussed above, it can give solutions 
of both energy and wavefunction of the ground state which may be the most important subject 
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to study in phenomenology of particle physics, and the solutions are almost-analytic, so that 
provide possibility for discussing physics of the concerned problems. Moreover, the simplic- 
ity of this method for dealing with complicated potential forms is remarkable. Therefore the 
Friedberg-Lee-Zhao method is absolutely applicable and advantageous. Of course it is worth 
further studies. 
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